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We study the weak-decay constants for the heavy pseudoscalar mesons, D, D s , B, and B s . For 
this purpose, we employ the extended nonlocal chiral-quark model (ExNLChQM), motivated by 
the heavy-quark effective field theory as well as the instanton-vacuum configuration. In addition to 
the heavy-quark symmetry and the nonlocal interactions between quarks and pseudoscalar mesons 
in ExNLChQM, a correction for the strange-quark content inside D s and B s is also taken into 
account and found to be crucial to reproduce the empirical values. From those numerical results, 
we obtain fD,D s ,B,B s = (207.53, 262.56, 208.13, 262.39) MeV, which are in good agreement with 
experimental data and other theoretical estimations. Using those numerical results, we estimate the 
CKM matrix elements and the Cabibbo angle with the various mesonic and leptonic heavy-meson 
decay channels, resulting in (|V cd |, \V CS \, \V u b\, \V td \/\V ts \) = (0.224, 0.968, < 5.395 x 10" 3 , 0.215) and 
0c — 12.36° which are well compatible with available data. 

PACS numbers: 12.15.Hh,12.39.Fe,12.39.Hg,13.20.Fc,13.20.He 
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I. INTRODUCTION 

Strongly interacting systems governed by quantum chromodynamics (QCD) manifests various interesting features. 
In terms of the quark and gluon degrees of freedom (d.o.f.), QCD is believed to indicate the confinement as well as 
the asymptotic freedom. At small strong-coupling regions, perturbative treatments have been successfully applied 
to investigate various physics problems at high energies. In contrast, due to the nontrivial QCD vacuum structure, 
the low-energy quark-gluon dynamics turns out to be very complicated, resulting in developments of lattice QCD 
(LQCD) simulations and effective approaches. Although LQCD has proved itself as a powerful method based on the 
first principle, i.e. QCD, it still contains several problematic issues, such as the notorious sign problem for finite- 
density QCD matters for instance. In turn, the effective approaches have been developed taking into account relevant 
symmetries in the nonperturbative regions and shown profound understandings for the nonperturbative QCD with 
many practical and successful applications 

Among the various kinds of those effective approaches, the instanton-vacuum configuration has accumulated reliable 
theoretical results and interpretations for numerous physical problems, such as the spontaneous breakdown of chiral 
symmetry (SBCS) [SHI], QCD vacuum properties [5UTT]. structure functions for the mesons and baryons [T2UT8] . QCD 
matter at finite temperature and/or density p~9j|20], an d so on. Note that the instanton is a semi-classical solution of 
the Yang-Mills equation in Euclidean space [2TJ [22] . The instanton vacuum is characterized by two phenomenological 
instanton parameters, i.e. average (anti) instanton size p « 1/3 fm and average inter- ( ant i) instanton distance R « 1 
fm [5HI]. In the dilute ensemble of these pseudoparticles, representing the nontrivial QCD vacuum, the quarks are 
delocalized and acquire their dynamically-generated effective masses, which depend on the transferred momentum [7]. 
It is worth mentioning that this momentum dependence plays the role of a UV regulator by construction so that 
artificial form factors for the quark-loop divergences are not necessary. 

By performing an appropriate bosonization process, we can obtain an effective chiral action (EChA), manifesting 
the nonlocal interactions between the quarks and pseudoscalar (PS) mesons [6]. It has been found that EChA is very 
useful to study the hadron physics in terms of the PS meson and quarks. There have been many attempts to extend 
EChA, defined properly in Euclidean space, to Minkowski space by considering an analytic continuation between 
them |23j|24]. For definiteness, we name the model defined in Minkowski space as the nonlocal chiral-quark model 
(NLChQM), although the naming is often used for different effective models in literatures. As shown in Ref. [23) [24]. 
NLChQM has applied to the studies for the light-cone PS-meson wave function and given reliable results in comparison 
with other theories and experiments. Note that NLChQM inherits the specific features of the instanton model as 
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understood, and possesses the similar symmetries and their breakdown patterns. 

Recently, much interest has been paid to the heavy meson and baryons, which contain the heavy-flavor quarks, i.e. 
c, 6, and t, from theoretical [25 -37 and experimental [38-45] points of view. Especially, energetic studies have been 
done by many experimental collaborations like D0 [38 , BaBar [39], Belle [40], BES [41], and so on. Also in the heavy- 
ion collision (RHIC) [42] and proton- (anti) proton scattering (LHC, Tevatron) [43H45] experiments, the heavy-hadron 
productions from the fire ball have been attracting interest. From the theoretical side, the nonperturbative QCD with 
the heavy-quark d.o.f. have been investigated actively in terms of the heavy-quark effective field theory (HQEFT) [46] . 
In HQEFT, there appear interesting features, such as the heavy-quark symmetry and velocity super-selection rule in 
the heavy-quark limit tuq —> oo. The heavy-quark symmetry consists of those for the heavy quark spin and flavor. 
According to them, we have the mass degeneracy between the pseudoscalar (scalar) and vector (axial-vector) heavy- 
meson states and the heavy-quark Dirac equation without the heavy-quark current mass. Successful applications of 
HQEFT can be found in Refs. [471453] and references therein. 

Considering the successes of NLChQM and HQEFT for the light- and heavy-quark sectors, respectively as mentioned 
above, it is quite natural and challenging to combine these two approaches and extend them to an effective model 
for the heavy-light quark systems. We want to name this new approach as the extended nonlocal chiral-quark model 
(ExNLChQM) [54 . Note that there were similar tries in the contexts of the Nambu-Jona-Lasinio model [55] and 
effective chiral Lagrangians [56j[57]. However, ExNLChQM inherits many useful and unique benefits from its origins, 
i.e. NLChQM and HQEFT, such as the heavy-quark symmetry, natural UV regulator, nonlocal quark-PS meson 
interactions, (relatively) strong constraints on the renormalization scale, and so on, as discussed previously. We 
found that ExNLChQM can reproduce the weak-decay constants for the non- strange heavy PS mesons for the SU(4) 
(u, d, c, 6)-flavor sector, fn and although there were uncertainties in the model parameters for the heavy-quark 
side [54]. It is worth mentioning that there are alternative attempts to incorporate the heavy-light quark systems, 
employing the instanton physics [58] [59] . 

Hence, in the present work, we want to extend our previous work into the strange heavy PS mesons for f]j s and 
fs s for the flavor SU(5) sector, i.e. (u,d,s,c,b) flavors. It is worth mentioning that the inclusion of the strangeness 
even for NLChQM in the flavor SU(3) is not a simple task, according to that the effects from the 1/N C corrections 
plays a crucial role. Thus, in order to incorporate the explicit strangeness with ExNLChQM, in principle, one needs 
to take into account the meson-loop corrections, corresponding to the 1/N C ones, as done in Refs. [60-62 . Moreover, 
it is worth mentioning that the axial-vector current conservation is of importance and the nonlocal contribution 
is necessary in order to conserve the current p~4j [24j [65] [66]. Instead of including the meson- loop and nonlocal 
contributions explicitly in the present work, however, we devise a very phenomenological way to circumvent these 
problematic issues for ExNLChQM, using the empirical information for « 132 MeV and fx ~ 160 MeV and their 
ratios. Detailed explanations will be given in Section III. 

Once having computed the theoretical results for /d,d s and /#,b s , we can estimate the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements via various leptonic and mesonic decay channels of the heavy PS mesons. Thus, these heavy 
PS-meson weak-decay constants are very important physical quantities for studying CP violation and have been 
studied extensively in numerous theoretical approaches, such as the light-cone formalism (LC) [25], light-front quark 
model (LQM) [26] , field-correlator method (FC) [27], Bethe-Salpeter method (BS) [28, 29 , relativistic quark model 
(RQM) [3D], QCD sum rule (QCDSR) [63 , and LQCD [3TH35] . From the numerical results in the present work, we 
obtain fD,D s ,B,B s = (207.53, 262.56, 208.13, 262.39) MeV, which are in good agreement with experimental and other 
theoretical values [25-37 . Using the present results, we estimate the CKM matrix elements and the Cabibbo angle, 
resulting in (|V cd |, I^UKbU^dl/l^l) = (0.224, 0.968, < 5.395 x 10" 3 , 0.215) and C = 12.36°. Again we find those 
numerical results are well compatible with available data [64]. These results indicate the reliability of the present 
theoretical framework. 

The present work is organized as follows: In Section II, we briefly introduce NLChQM motivated by the instanton 
vacuum configuration. NLChQM is extended to ExNLChQM with HQEFT in Section III. Section IV is devoted to 
numerical results and corresponding discussions. Summary of the present work and future perspectives are given in 
Section V. 



II. NONLOCAL CHIRAL-QUARK MODEL (NLCHQM) 

First, we want to make a brief explanation on an effective model, motivated by the instanton vacuum configura- 
tions. The model describes the highly-nonlocal interactions between the quarks and nonperturbative gluons, i.e. the 
instanton solution in Euclidean space [5HZ]. By a proper bosonization process similar to that for the NJL model [U[2], 
integrating out all the meson fields except for the SU(7V/) multiplet PS mesons, one arrives at the following effective 
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FIG. 1: (Color online) Effective light-quark mass M q as a function of the momentum transfer p, using Eq. |2]) (solid) and 
Eq. Q (dot). Here, we employ M g , = 350 MeV and A q = l/p = 600 MeV. 



action: 

S eS [(j),m q ] = -Spin ifi- 



■iy/W q Usy/M q 



(1) 



where <j) and rh q stand for the light PS-meson field and current-quark mass matrix for the flavor SU(3), 
diag(m n , m^, m s ). Sp and M q = M q (d 2 ) denote the functional trace over relevant spin indices and momentum- 
dependent effective-quark mass for the light-flavor quarks q = d, s), respectively. Generically, M q is defined by the 
modified Bessel functions in momentum space by [7] 



M q (p 2 ) = M qfi F 2 (p 2 ) 



2t 



/ (i)i^i(i) - h(t)K (t) - -him^t) 



p 2 p 



(2) 



We note that F(p 2 ) comes from the Fourier transform of the instanton zero mode [7]. In many applications, M q is 
parameterized for simplicity as follows: 



M q {p 2 ) = M qfi 



2A 2 

q 



2A2+p2 



(3) 



In the present work, we will use Eq. ([3| for numerical calculations, since there are only small differences between the 
usages of Eqs. ^ and ^ as shown in Figure]!] The self-consistent equation of the model for the chiral limit reads [7_: 
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d 4 p M 2 q 



) 4 p 2 + M 2 ' 



(4) 



If we make use of the phenomenological instanton parameters R ~ 1 fm and p ~ 1/3 fm, the value of M q ^ is 
determined to be about 350 MeV via Eq. Thus, we will use this value for M q ^ and a renormalization scale of the 
model A q « 1/p « 600 MeV throughout the present work. 

Taking into account the nonzero current-quark mass beyond the chiral limit, the weak-decay constant for the PS 
meson, we have the following formula for the leading local contribution for f n and fx- 



4N r 



d 4 p N f N g ^MjW g M f + M g - f(M> + M r g ) 



(p 2 +M / 2 )( p 2 + ^2) 



(5) 



where we have defined the notations as M q = M q + m q and M' = dM/d\p\. The subscripts / and g denote the quark 
for each flavor inside the PS meson <j) ~ fg. Throughout this work, we employ the normalization as f n = y/2F^^ where 
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F n « 93 MeV. Note that we introduced a phenomenological multiplicable factor M q and will explain the meaning of 
this parameter in detail below. First, choosing A// = M g = 1, we have the following numerical values using Eq. ([§: 

113 MeV for m uAs = (5, 5, 100) MeV. (6) 

In order to obtain the empirical value for f n « 132 MeV, one has Af Uj d = 1.17, taking into account the isospin 
symmetry m u = rrid as in Eq. As for the kaon, J\f s = 1.71 reproduces the empirical value for fx ~ 160 MeV. 
Below, we want to explain the physical meanings of these Aff parameters in detail: 

• In order for satisfying the axial- vector- current conservation for the PC AC relation in the nonlocal effective model 
in terms of the flavor SU(3) symmetry, one needs to include nonlocal contributions additionally to Eq. ^ p~4j 
[24] [65] [66] , although we have ignored them for simplicity. They provide about (10 ~ 30) % increase of relevant 
physical quantities in this model. Hence, N u ,d slightly bigger than unity compensate the absence of the nonlocal 
contributions in Eq. This simplification can be understood easily by the following: 

n = f*+fT j ^M a Mp%, (7) 

where the superscripts L and NL indicate the leading local and sub-leading nonlocal contributions, respectively, 
while the subscripts a and j3 represent the quark flavors inside the PS meson in the flavor SU(3) symmetry 
m u ~ m d ~ m s <C A q . 

• If we go beyond the flavor SU(3) symmetry and the heavier strange quark comes into play, breaking the symmetry 
(mu,d <C ra s ), the situation gets complicated: One needs to consider the 1/N C corrections [I0l|62], which also 
corresponds to the meson- loop corrections (MLC). Therefore, A/" s , which is considerably larger than unity, plays 
the role of MLC effectively as 

U. « ATcMpfi. + /r° ^4. (° r W)' ( 8 ) 

Here, the subscript <p s stands for the PS meson including the heavier strange quark, i.e. the kaon. 

Thus, M q can be understood as a very phenomenological compensation for the above two simplifications to reproduce 
the empirical data. The differences in N u ,s,d m &y also indicate the distinctive interaction strengths between the 
flavors and the QCD vacuum, as observed in usual constituent-quark models. In what follows, N u , s ,d will be used 
for computing the heavy- meson decay constants with strangeness, instead of cosidering the nonlocal (F^) and MLC 
(F| /[LC ) contributions explicitly. We verified that this phenomenological approach makes the problems much simpler 
to a great extent analytically as well as numerically. Moreover, the numerical results obtained from this approach are 
turned out to be qualitatively-well compatible with other experimental and theoretical studies as shown in Section 
IV. 
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m s 
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M D 


M Ds 


M B 
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5 


5 


100 


1270 


4670 


1867 


1968 


5279 


5366 



TABLE I: Numerical inputs for the relevant masses [MeV] for the calculations [64] . 



III. EXTENDED NONLOCAL CHIRAL QUARK MODEL 

Now, we are in a position to introduce ExNLChQM as derived in the previous work [54] . We note that this model 
is in principle equivalent to that suggested in Ref. [55] in many aspects, besides several specific features, i.e. quark-PS 
meson nonlocal interactions and natural UV regulator by construction for instance. Those features are inherited from 
the instanton- vacuum effects. By taking into account the (u,d,s,c,b) flavors and the structure of Eq. 0, one can 
construct the following effective chiral Lagrangian for the heavy (Q) and light (q) quark system in Minkowski space: 

£ ExNLCh Q M^ m] = ^ _rn q -m Q - VM^UsVm] iP, (9) 

where stands for the heavy (H) and light (L) mesons. rh q ^Q indicate the current-quark mass matrices for the light- 
and heavy-flavor quarks for q = (u,d,s) and Q = (c, b): rn q = diag(ra n , m^, ra s , 0, 0) and rh q = diag[0, 0, 0, m c , m^). 
The quark spinor is assigned as ip = (u, d, s, c, b) T . Although the effective Lagrangian in Eq. (I9J) has a SU(5) 
symmetric form for the flavors, the symmetry is broken explicitly by the difference between the current-quark masses 
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as in Table [T| Moreover, the flavor SU(5) symmetry is broken further by the distinctive effective quark masses for the 
heavy and light quarks as will be discussed later. Note that the chiral symmetry is explicitly broken as understood by 
Eq. ri9 and SBCS is emerged by finite M values, corresponding to the effective quark mass, A4 = diag(M g , Mq) = 
diag(M n , Md, M s , M c , Mb) as in Eq. Q [54]. We define the effective heavy-quark mass, assuming it possesses the 
same analytic structure with that for the light quark, in Minkowski space: 



M Q (p 2 ) =Mg, 



2A 2 Q 



2Aj 



(10) 



where Mq^ indicates the effective heavy-quark mass at zero virtuality. Aq denotes a renormalization scale for the 
heavy quarks. The nonlinear heavy-light PS-meson fields for the SU(5) symmetry can be constructed as follows: 



Us = exp 



(11) 



Here, the weak-decay constant is assigned by the normalization F<$> = as mentioned. The explicit form of the 

24-plet U5 matrix is given in Appendix. This sort of the extension of a flavor group is in principle equivalent to that 
suggested in Ref. [67] . 

Since we are interested in computing the weak-decay constants in the present work, we expand the nonlinear PS- 
meson fields in the effective Lagrangian up to for the current oc (0| Jw\®)- Then, the effective Lagrangian in 
Eq. (9| can be represented in three separate parts, i.e. light-light (LL), heavy- heavy (HH), and heavy-light (HL, LH) 
quark terms as done in Ref. [54] : 



^ExNLChQM _ r LL , /-HH 



q(ip q 

1 



m„ 



r (HL,LH) 
~eff 

M q )q- 



-qy/M q [i 75 L] y/M q q 



LL 



[0 (ip Q -m Q -M Q )Q} 



HH 



F H 



Q^Mq [i lb H] ^M q q 



J HL 



F H 



qy^Wq [iHj 5 ] ^/MqQ 



(12) 



LH 



where L and H denote the light and heavy PS-meson fields defined in Eq. (32) in Appendix. Note that, however, we 
have not considered the QQ-meson states, corresponding to the quarkonia states, in the above effective Lagrangian. 
In Eq. (12), we have also taken into account the effective heavy-quark mass term, Mq in £™, in addition to the 
current mass tuq. According to HQEFT, one can redefine the heavy-quark and heavy-meson fields as follows [55] : 



Q(x) 



-imav-x 



Qv(x), H = e 



-imov-x 



H v , H = e 



imov-x 



H v . 



(13) 



Note that v is the heavy-quark velocity. Here, we choose v — (+1,0,0,0) for definiteness, representing the heavy 
quark at rest [46 . It is an easy task to redefine C™ in Eq. (12) using the Dirac equation for the heavy-quark field 
in Eq. (13), resulting in 



/-HH 



] v [f{iv • d) - M Q ] < 



)v [(iv ■ d) - Mq] . 



In the second step of Eq. (14), we use the velocity projection for the heavy quark, fQ v 
Eq. (13), we can rewrite £^^' LH ^ in the following form: 



(14) 

Q v [46]. Similarly, using 



/-HL 



-eff 



F H 



Qv\[Mq 



i + i> 



/-LH 

^eff 



F H 



'M„ 



(iH v ) 



-75 



(15) 



Making use of a generic functional integral technique for the Grassmann variables given for the two Grassmann 
variables q and Q v [54], finally, we can arrive at an EChA for the heavy-light quark systems from the effective 
Lagrangian density, given in Eq. (12): 



C LL+HL+LH 
^eff 



— iSp In 



ip-M q 



1 

F~l 



F H 



M Q H^yW q )(ivd-M Q )- 



1 

Fh 



'M q R^ 



'Mr 



.(16) 



It is worth noting that the effective action in Eq. (16) is in principle equivalent to the first term of Eq. (36) in Ref. [55] , 
except for the momentum dependent quark-PS meson coupling strengths. 
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As a next step, we estimate Mqq in Eq. (10) from a simple phenomenological analysis. In this consideration, the 
heavy PS-meson mass can be understood as 



M, 



H 



K + M,, ] L + [mg + M, 



Q,oJ H 



(17) 



where we have ignored the binding energy for the meson. From the experimental data for D and B mesons, we can 
write 



M B 



1869.57 MeV 
5279.17 MeV 



(ra c 



m q 



M q , 



M Q,o) 



1625.0 MeV 
5025.0 MeV 



Mg, M Q)0 
M Q)0 Mq, 



244.57 MeV 
254.17 MeV. 



(18) 



The numerical inputs, which are taken from Ref. [64 , are summarized in Table [TJ As shown in Eq. (18), it is necessary 
to add the effective heavy-quark mass to reproduce the heavy-meson mass appropriately. 

Here are some explanations on this additional heavy-quark mass: In the instanton model for example, the light 
current quarks obtain their momentum-dependent effective masses via the nontrivial interactions with the instanton 
ensemble, which represents the nonperturbative QCD vacuum [5 -7 . This mechanism is also equivalent to SBCS, 
resulting in the nonzero values for the chiral condensate (qq) and various low-energy constants, such as F$. Hence, 
it is reasonable to assume a similar mechanism for the additional heavy-quark mass Mq ? o in Eq. (18), i.e. Mq$ 
is considered to be originated from the highly nontrivial interactions between the heavy quark with tuq and the 
nonperturbative QCD vacuum. To estimate Mq ? q, we have chosen M Qj q « 350 MeV as in Section II. If we consider 
the binding energy for the mesons, the estimated value Mq ; o = (240 ~ 250) MeV must be the lower bound of 
its real one. From these observations and previous discussions, we summarize our educated assumptions from a 
phenomenological point of view: 

• The finite effective heavy-quark mass, Mq^ is generated from the similar mechanism with that for the light 
quark-instanton interaction, representing the nontrivial QCD vacuum effects. For the numerical calculations, 
Mq j0 can be estimated by Eq. (17). 



• If the instanton ensemble is not affected much by the heavy sources such as the heavy quarks, one uses Aq « 600 
MeV again for Eq. ( 10 ) like that for the light quarks. However, we will see that the value of Aq should be changed 
to reproduce the experimental data for appropriately. The difference between A q and Aq may signal the 
broken flavor SU(5) symmetry. 



Thus, we will use the effective heavy-quark mass in the momentum space as in Eq. (10), similar to the light-quark 
case. 

The heavy- meson weak-decay constants is defined as follows [55] : 



<0|g(a;)7„(l - j 5 )Q v (x)\H(p)) = ipJ H , 



(19) 



where p stands for the on-mass shell momentum of H with the velocity v. In order to evaluate the matrix element in 
Eq. (19), we employ the external-field method as done in Ref. [54]. After a straightforward functional manipulation, 
we obtain the following formula for the heavy- meson weak-decay constant: 



2N c M q f k 2 dk ^/M qfi M Qfi (AM H - M Qfi )[k 2 + (M o , - AM H ) 2 + 2A 2 } 4 
AM H J 7T 2 [k 2 + (M Qfi -AM H ) 2 + 2A 2 g ]*[k 2 + (M Qfi -AM H ) 2 ] + (2A 2 g )*M 2 



2A 2 



[k 2 + (M Qi0 - AM H ) 2 + 2A2] 



2A| 



(k 2 + 2A 2 ) 



(20) 



where, from the on-mass shell condition for the meson H at v [55], one is lead to v • p ~ AMh = Mh — tuq. 



IV. NUMERICAL RESULTS 



In this Section, we present the numerical results and related discussions. First, we show them for /d,d s (left) and 
Ib,b s (right) for M q = 1 in Eq. (20) in Figure [2] as functions of Aq at A q = 600 MeV. The curves for the strange 
and non-strange mesons are drawn in the solid and dot lines, respectively. Horizontal lines denote the experimental 
data, while the errors are represented by the shaded areas. The experimental data are taken from Refs. [36j[37]. As 
shown in Figure [2j for the non-strange mesons, i.e. /£>,#, we see that the data are well reproduced when Aq « 1 
GeV: /d,b = (207.54, 208.13) MeV, where the experimental data provide f DyB = (206 ± 8.9 204 ± 31) MeV [311 [37]. 
Although one may choose J\f u ,d — 1-17 as done for f n in Section II and modify the Aq value from 1 GeV, we assume 
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FIG. 2: (Color online) f D , Ds (left) and f B ,B s (right) as functions of Aq at A g = 600 MeV for Af u , d ,s = 1 in Eq. The 
numerical results for the non- strange and strange mesons are shown with the solid and dot lines, respectively. The experimental 
data are taken from Refs. 36, 37 . Note that the experimental errors are given with the shaded areas. 



that the effects from 5Af u ,d, which amounts the deviation from unity, i.e. J\f u ,d = 1 H-AA/^d, are absorbed into Aq for 
brevity. 

Physically, this larger renormalization scale for the heavy quark Aq « 1 GeV can be understood intuitively as 
follows: As for the light meson, consisting of the two light quarks, the typical scale has been determined in general by 
A q « 600 MeV as in many instanton-motivated models like the present work. However, taking into account the heavier 
degree of freedom for the heavy mesons, including the heavy quarks tuq ~ a few GeV, one must introduce larger scale 
in comparison to the light mesons. Hence, this effect can be considered effectively as well as phenomenologically in 
the larger renormalization scale to reproduce the data. 

Taking this strategy, we fix Aq = 1 GeV for every heavy-quark flavors from now on. The numerical results for 
fn 3 turn out to be much smaller than the experimental value, fo s = 200.78 MeV with Aq = 1 GeV, whereas 
we do not have experimental data for fs s - These underestimations for the strange mesons can be understood by 
the same problem for fx with the explicit flavor SU(3) symmetry breaking as discussed in Section II. Hence, if we 
chose Af s = 1.71, which was introduced to reproduce the empirical value for the numerical results for /d s ,b s are 
drastically improved as shown in Figure [31 Numerically, we have /d s ,b s = (262.56, 262.39) MeV to be compared with 
the experimental data fu s — (257.3 ± 5.3) MeV. From these observations, the correction factor J\f s for the strangeness 
in the present model is crucial, and we have a tendency approximately from the present model: 

(/D~/s«210MeV) < (f Da w/ B . « 260 MeV). (21) 

Our final numerical results for /h,h s and their ratio between them are listed in Table |ll| 





Id 


Sb 


fo s 


f Bs 


fnjfn 


SbJSb 


Present 


207.54 


208.13 


262.56 


262.39 


1.27 


1.26 


Experiment 


206 ±8.9 


204 ± 31 


257.3 ± 5.3 




1.27 





TABLE II: Numerical results for f D ,D s ,B,B s [MeV] and their ratios, computed for A q , Q = (600, 1000) MeV, M q , = 350 MeV, 
and J\f u ,d,s = (1, 1, 1.71). The experimental data are taken from Refs. [36l [37] , 

In Figure [4j we compare our numerical results with other theoretical estimations for fn,D s (left) and /b,b s (right). 
Other theoretical results are taken from the clover- improved quenched LQCD [31], field-correlator method (FC) [27] . 
QCD sum rule (QCDSR) [63], light-front quark model (LQM) [26], Bethe-Salpeter method (BS) [28] [29], relativistic 
quark model (RQM) [30], and light-cone wave function (LC) [25]. As for /o, all the theoretical results are quali- 
tatively compatible with the data, while that from BS [28] [29] overestimates considerably. In contrast, LQCD [3T] 
underestimates the data for /^ s , depending on the simulation schemes. As shown in the right panel of Figure [4] the 
present result for fs are in good agreement with the data which contain huge error bars. All the theoretical results 
are still comparable with the experimental center value for fg = 204. The situation changes, however, drastically for 
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FIG. 3: (Color online) fD,D s (left) and fB,B s (right) as functions of Aq at A q = 600 MeV for Af u ,d = 1 and Af s = 1.71 in 
Eq. (20). The numerical results for the non-strange and strange mesons are shown with the solid and dot lines, respectively. 
The experimental data are taken from Refs. [36, 37 . Note that the experimental errors are given with the shaded areas. 



Jb s - Note that the present result and that from LC [25] give Jb s > 260 MeV, while others locate below about 240 
MeV. Hence, the experimental measurements for fs s can be a good place to test theoretical reliability of each model 



as shown here. Numerical values in comparison with other theoretical calculations are summarized in Table III 





Present 


LC [25] 


LQM [26] 


FC [27] 


BS EH [29] 


RQM [30] 


LQCD [31] 


LQCD [32] 


LQCD [32] 


LQCD [34] 


LQCD [35] 


Id 


207.53 


(206 ± 8.9) 


211 


210 ± 10 


230 ± 25 


234 


211±? 2 ± 14 


194±^ ± 10 


218.9 ± 11.3 


192 ± 30 


197 ±9 


fn s 


262.56 


267.4 ± 17.9 


248 


260 ± 10 


248 ± 27 


268 


23llo ± 12 


2131^ ± 11 


260.1 ± 10.8 


214 ±33 


244 ±8 


fB 


208.13 


(204 ±31) 


189 


182 ±8 


196 ± 29 


189 


179±| 6 ± 18 


164±{l ± 8 


196.9 ±8.9 


171 ± 29 






262.39 


281 ± 54 


234 


216 ±8 


216 ±32 


218 


204±g 8 ± 16 


185^ 3 ± 9 


242 ± 9.5 


193 ± 32 





TABLE III: Theoretical results for fD,D s ,B,B s from the present calculations, light-cone wave function (LC) [25 , light-front 
quark model (LQM) [26], field- correlator method (FC) [27 , Bethe-Salpeter method (BS) [28] [29], relativistic quark model 
(RQM) [30], clover-improved quenched LQCD (31 j , 0(a) -improved light-quark action [32], (2± l)-flavor asqdat action [32], full 
one-loop correlator [34], and Symanzik gauge action [35] . 



In what follows, we extend our discussions to very rough estimations for the heavy- meson nonrelativistic wave 
function and CKM matrix elements, based on the numerical results obtained. In the nonrelativistic quark models, 
the heavy PS-meson decay constant is given by the following formula [68] : 
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w 



H 



l*tf(0)| 2 , 



(22) 



where \I/#(0) is the wave function for the heavy meson at zero separation between the heavy and light quarks inside. 
Using the computed values for fjj as shown in Table [Til we obtain 



\^ D,D S ,B,B S (0)| = (0.08,0.11,0.14,0.18), 
which result in the following ratios, 



hMQ)l 
1*^.(0)1 



= 0.73, 



I*b(0)| 
I*b.(0)| 



= 0.78,, 



I*d(0)I 
I*b (0)| 



= 0.57, 



\*dM\ 

l*B.(0)| 



= 0.61. 



(23) 



(24) 



Note that the ratio between D(D S ) and B(B S ) is considerably smaller than unity, being different from the rough 
assumption made in Ref. [69], whereas the estimations for the ratio between D(B) and D S (B S ) are relatively similar 
to each other. These numerical results for the ratios will be informative to construct and estimate the heavy- meson 
wave functions in the nonrelativistic quark models. 
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FIG. 4: Theoretical results for fD,D s ,B,B s from the present calculations, clover-improved quenched LQCD 31 , field- correlator 
method (FC) [27], QCD sum rule ggg,' light-front quark model (LQM) [26], Bethe-Salpeter method (BS) [28] [29], relativistic 
quark model (RQM) [30 , and light-cone wave function (LC) [25] . 



Now, we want to apply the numerical results for fD,D s ,B,B s to some phenomenological problems. As known well, 
the PS-meson leptonic decay is a good subject to study the weak interactions corresponding to the CKM matrix 
elements. Here, fjj plays a critical role. The partial decay width for P — ^ Iv reads: 

T{P^tv) = ^f P M P M 2 (l-g)V /s | 2 , (25) 

where Gp, P, /p, Mp, and Mi denote the Fermi constant, generic PS meson, PS-meson weak-decay constant, PS- 
meson mass, and lepton mass, respectively. Vf g stands for the CKM matrix element for the quark flavors / and g. 
As for D + — )> from the CLEO experiment, the branching ratio B(D+ — >• ii^v) was given by (3.82 ± 0.32 ± 

0.09) x 10 -4 [64]. Note that £> is defined by TpT(P — >• iv)/h, in which Tp denotes the life time of the PS meson. 
Using the numerical result for /p> given in Table \H\ we have the CKM matrix element \V c d\ = 0.224, whereas it 
becomes 0.225 for the experimental fp, value. Similarly, as for Df and P + from the present (experiment) results, 
we have \V CS \ = 0.968(0.975) and \V ub \ < 5.395(5.501) x 10- 3 (90%CL), using S(P>+ -> p+v) = 5.8 x 10~ 3 and 
B(B+ fjL+v) < 1.0 x 10~ 6 [64]. Note that the PDG fits provide \V cd \ = 0.2252 ± 0.0007, \V CS \ = 0.97345±Hooi5> 
and \V u b\ = (3.47±g;J|) x 10~ 3 . As understood, our numerical results are in good agreement with the PDF fits. All 
the numerics for the CKM matrix elements are listed in Table [El 

The mesonic decay of B° can be computed by the following formula [69 : 



T(B° -> U+P-) = ^-\V cb \ 2 \V fg \ 2 M% f% \Hwl)\ 2 \J Hh ^ B) {{l ~ VA)[(l + VA)~ B]}\ (26) 

where we have used the notations A = M^/Mp and B = M|/M|, and \(x,y,z) indicates the following kinematic 
function, i.e. 

X(x, y, z) = [x 2 -(y + zf] [x 2 - (y - z) 2 }. (27) 
The universal form factor Tp for the PS meson can be parameterized as follows: 



a s(m 2 b ) 
a s (m 2 c ) 



1 



1 - w 2 p /wl 



2 • 



(28) 



Here are two dimensionless quantities: wp = [Mp — (Mp — Mp>) 2 )/ (MpMp>) and wq corresponds to the cutoff mass 
for the form factor parameterization which can be determined by experiments and/or theories. In Ref. [69], wo 
was estimated as 1.12 ±0.17 that we will use for further discussions below. We also choose a s (m 2 ) = 0.189 and 



a s (m 2 ) = 0.29 [69]. Using Eq. (26) for P = (D ,D S ), we can write a ratio as follows: 



B(B° ^D+D ) _ ^ x tan2 Qc x /p_ \ 2 = Q 7g x tan2 ^ (29) 
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where tan 6 c = \V us \/\V u d\ with the Cabibbo angle 6c- Considering the experimental data for the branching ratios 
without mixings, B(B° D + D~) = (2.11 ±0.31) x 10~ 4 and 23(5° -> D+D~) = (7.2 ±0.8) x 10" 3 [64], one obtains 
6 C ~ 11.11°. Note that, if we choose F{w 2 D ) « ^(wj^J in Eq. p6| as in Ref. [70 , # c becomes about 12.36°. These 
estimations are compatible with the empirical value 6c ~ 13° |64j . 

The CKM matrix elements with t flavor can be indirectly obtained by the following equation [64] : 



SB 



ts 



v td 



(30) 



where Am# 5 £ s stand for the mass difference between the heavy and light mass eigenstates of the B°-B mixing in 
the standard model (SM), so called as the oscillation frequencies [64], in which those values are given as Ara# = 
(0.507 ± 0.005) ps _1 and Am Bs = [17.77 ± 0.10 (stat) ± 0.07 (syst)] ps -1 . £ B is the flavor SU(3) symmetry breaking 
factor, defined by 



b b s fB a 
Bb Jb 



(31) 



Here, Bb,b s denote the bag parameter appearing in the SM calculations for the B° — > B° transition. If we set their 
ratios to be unity as done in Ref. [70], we have from the present numerical results = 1.26, which is well compatible 
with that given in the LQCD simulation [71], £b,lqcd = 1.2lt°;°^. Ref. [72] gives £ B = 1.243 ± 0.021 ± 0.021 via 
a global fit using the LQCD data. Plugging the empirical values for Am B ^ Bs and Mb,b s into Eq. (30) and using 
the numerical result for we have |Vtd/Vt s | = 0.215, which is similar to the PDG fit, ~ 0.214. Note that a LQCD 
simulation [72] provides |Vtd/Vt s | = 0.211 ±0.001 ± 0.005 [64 . The numerical results are summarized in Table IV 





\v cd \ 


\v cs \ 


\V ub \ x 10- 3 


\Vtd/Vts\ 


Present 


0.224 


0.968 


< 5.395 


0.215 


PDG fit 


0.2252 ± 0.0007 




o 47 +o.i6 

° -0.12 


0.211 ± 0.001 ± 0.005 (LQCD [61 [72]) 



TABLE IV: CKM matrix elements including the heavy-quark flavors, \V c d\, \V CS \, and \V u b\ 



V. SUMMARY AND PERSPECTIVES 



In the present work, we have investigated the weak-decay constants for the heavy PS-mesons consisting of the 
(iz, s,c, b) flavors, i.e. fD,D s ,B,B s - To this end, we employ flavor SU(5) ExNLChQM which was developed being 
motivated by HQEFT and NLChQM. We chose the cutoff masses for the light and heavy quarks differently as A q = 600 
MeV and Aq = 1000 MeV, respectively, considering that the nontrivial QCD vacuum effects are finite even for the 
heavy quark inside the meson. In what follows, we list important observations in the present work: 

• We employ a phenomenological correction factor for the inclusion of the strange quark, compensating the 1/N C 
corrections. By doing that, we obtain fD,D s ,B,B s — (207.53, 262.56, 208.13, 262.39) MeV, which are qualitatively 
compatible with available experimental and theoretical values. Thus, this correction factor turns out to be crucial 
to reproduce the data. 

• In other words, this qualitatively good agreements with the data indicate considerable contributions for the 
heavy quarks from the nontrivial QCD-vacuum effects. Moreover, it is justified that the vacuum structure due 
to the heavy quarks is modified in comparison to that for the light quarks by seeing that A q ^ Aq, since A 
value denotes the average size of (anti)instantons. This tendency also signals the breakdown of the SU(5) flavor 
symmetry at the quark level, in addition to m q ^ tuq. 

• The nonrelativistic heavy- meson wave functions at zero separation between the light and heavy quarks, \I/#(0), 
are calculated. From those numerical results, we observe that the ratio between the wave functions for the strange 
and non-strange heavy mesons is less than unity. This observation is quite different from simple consideration 
in quark models, i.e. |^#(0)|/|^# s (0)| = 1. 

• The CKM matrix elements are estimated by using the various leptonic decays of the heavy PS mesons, resulting 
in that (|V^|, \V C3 \, \V u b\) = (0.224,0.968, < 5.395 x 10 -3 ). Taking into account the mass difference between the 
heavy and light mixing of B°-B° states, although there are uncertainties arising from the bag parameters in SM, 
we have |V^|/|^s| = 0.215. Again, these values are well comparable with their empirical and LQCD simulation 
data. 
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• Accounting for the hadronic decays for the heavy meson and their ratios, we compute the flavor SU(3) symmetry 
breaking factor as = 1.26. The Cabibbo angle can be also estimated using this information, resulting in 
9c = 12.36°, depending on the choice of the PS-meson form factor. We note that these values are consistent 
with its empirical data and LQCD results. 

Taking into account all the results summarized above, we can conclude that the present theoretical framework, 
ExNLChQM is a qualitatively successful model to study the heavy PS meson, which consist of the heavy and light 
quarks. As for the future improvements of the present form of ExNLChQM, 1) the nonlocal contributions for con- 
serving the axial- vector current can be included to compute the weak-decay constants. Furthermore, 2) the 1/N C 
corrections can be included in the present framework in principle via standard functional treatments, instead of using 
the phenomenological correction factor, although this task must quite challenging. Nonetheless for these possible 
improvements, from a phenomenological point of view, the present theoretical framework is still a considerably use- 
ful tool to investigate various physical quantities for the heavy-light quark systems, such as the Isgur-Wise function, 
heavy-light meson coupling constants, heavy-meson effective Lagrangians, and so on. Related works are under progress 
and appear elsewhere. 
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Appendix 



The SU(5) PS-meson fields are defined as follows: 



U 5 = v / 2 
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